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Using Taylor expansions of the pressure obtained previously in studies of 2-flavor QCD at non¬ 
zero chemical potential we calculate expansion coefficients for the energy and entropy densities up 
to Cl(Atq) in the quark chemical potential. We use these series in p,q/T to determine fines of constant 
entropy per baryon number {S/Nb) that characterize the expansion of dense matter created in 
heavy ion collisions. In the high temperature regime these lines are found to be well approximated 
by fines of constant piq/T. In the low temperature phase, however, the quark chemical potential 
is found to increase with decreasing temperature. This is in accordance with resonance gas model 
calculations. Along the fines of constant S/Nb we calculate the energy density and pressure. Within 
the accuracy of our present analysis we find that the ratio p/e for T > To as well as the softest point 
of the equation of state, {p/e)min — 0.075, show no significant dependence on S/Nb- 

PACS numbers: 11.15.Ha, ll.lO.Wx, 12.38Gc, 12.38.Mh 


I. INTRODUCTION 

Recently studies of QCD thermodynamics on the lattice have successfully been extended to non-vanishing quark 
(or baryon) chemical potential using Taylor expansions nil as well as reweighting techniques Q around the limit 
of vanishing quark chemical potential {ptq = 0) and analytic continuations of numerical calculations performed with 
an imaginary chemical potential 0. The complementary amroach at fixed baryon number mi has also been used 
recently in simulations of 4-flavor QCD with light quarks [ 3 . These calculations yield the additional contribution 
to the pressure that arises from the presence of a net excess of baryons over anti-baryons in a strongly interacting 
medium. We will further explore here the approach based on a Taylor expansion of the pressure and make use of our 
earlier calculations for 2-flavor QCD H to determine also the energy and entropy densities at non-vanishing baryon 
chemical potential. 

In a heavy ion collision a dense medium is created which after thermalization is expected to expand without further 
generation of entropy (S') and with fixed baryon number (Nb) or, equivalently, with fixed quark number Nq = SNb- 
During the isentropic expansion the ratio S/Nb thus remains constant^. The cooling of the expanding system then is 
controlled by the equation of state on lines of constant S/Nb- From a knowledge of the energy density and pressure 
at non-vanishing quark chemical potential we can calculate trajectories in the /iq-T phase diagram of QCD that 
correspond to constant S/Nb and can determine the isentropic equation of state on these trajectories. 

This paper is organized as follows. In the next section we summarize the basic setup for calculating Taylor 
expansions for bulk thermodynamic observables and present a calculation of the Taylor expansion coefficients for 
energy and entropy densities at non-vanishing quark chemical potential up to 0{pq)- We use these results in Section 
III to determine lines of constant S/Nb and study the temperature dependence of pressure and energy density along 
these lines. Our conclusions are given in section IV. 


II. TAYLOR EXPANSION OF PRESSURE, ENERGY AND ENTROPY DENSITY 

The analysis we are going to present in this paper is based on numerical results previously obtained in simulations of 
2-flavor QCD on lattices of size 16^ x 4 with an improved staggered fermion action HQ In these calculations Taylor 
expansion coefficients for the pressure have been obtained up to 0{pP) for a fixed bare quark mass value (m = 0.1) 
which at temperatures close to the transition temperature (Tq) corresponds to a still quite large pion mass of about 


^ Here Nb denotes the net baryon number, i.e. the excess of baryons over anti-baryons. 
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770 MeV. However, in particular at temperatures above the QCD transition temperature the remaining quark mass 
dependence of thermodynamic observables is nonetheless small as deviations from the massless limit are controlled 
by the quark mass in units of the temperature, which is a small number. 

We closely follow here the approach and notation used in Ref. . We start with a Taylor expansion for the pressure 
in 2-flavor QCD for non-vanishing quark chemical potential (and vanishing isospin chemical potential, /r/ = 0), 

- OO 

^ , ( 1 ) 

n—0 


with 


c„(r, mq) 


1 1 5"lnZ(T,/x,) 

'W\VT^ d{^iq/TY 


( 2 ) 


Here we explicitely indicated that the Taylor expansion coefficients c„ depend on temperature as well as the quark 
mass^. In the context of lattice calculations it is more customary to think of this quark mass dependence in terms 
of a dependence on an appropriately chosen ratio of hadron masses characterizing lines of constant physics. We thus 
may replace niq by a ratio of pseudo-scalar (pion) and vector (rho) meson masses, mq = (mps/mvY. 

Due to the invariance of the partition function under exchange of particle and anti-particle the Taylor expansion is a 
series in even powers of Hq/T] expansion coefficients c„ vanish for odd values of n. From the pressure we immediately 
obtain the quark number density. 


nq 
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Using standard thermodynamic relations we also can calculate the difference between energy density (e) and three 
times the pressure, 


€ — 3p 

2^4 
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where 


c'n{T,mq) = T 


dc„(r, mq) 


^ dT 

Combining Eqs.mS] and^we then obtain Taylor expansions for the energy and entropy densities. 
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( 6 ) 


The expansion coefficients Cn{T,mq) have been calculated in 2-flavor QCD at several values of the temperature and 
for a fixed value of the bare quark mass, m = 0.1 Q. We note that the bare coupling th introduces an implicit 
temperature dependence if m is kept fixed while the lattice cut-off is varied. The latter controls the temperature of 
the system on lattices of temporal extent Nt, i.e. T~^ = NtO. The difference between derivatives taken at fixed th 
and fixed m^, however, is negligible at high temperature, where the quark mass dependence of bulk thermodynamic 
quantities is 0{{mq/T)‘^) and, of course, vanishes at all temperatures in the chiral limit as this also defines a line of 
constant physics. We thus approximate the derivative, Eq. by a derivative taken at fixed rh which becomes exact 
in the chiral limit. These derivatives are evaluated using finite difference approximants. Alternatively we may express 
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We assume that the thermodynamic limit has been taken so that we can ignore any volume dependence of thermodynamic quantities. 
For an aspect ratio = 4 and the quark mass value used in our calculation finite volume effects are expected to be small. 
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derivatives of Taylor expansion coefficients with respect to T in terms of derivatives with respect to the bare lattice 
couplings (3 = Qlg^ and m and rewrite Eq. I^as, 


d/3 dcn{T, rriq) _ dm dCn{T, mg) 
da dp da drh 


Here the two lattice /3-functions, ad/3/da and adm/da, have to be determined by demanding that the temperature 
derivative is taken along lines of constant physics |l(i lll| . The /3-function controlling the variation of the bare quark 
mass with the lattice cut-off, adm/da, is proportional to the bare quark mass m and thus vanishes in the chiral limit 
da. For small quark masses the first term in Eq. [3 thus gives the dominant contribution to c^(T, m^). We have 
evaluated this first term in Eq. [3 for n = 2 and find that this approximation for c '2 agrees within errors with the 
calculation of c '2 from finite difference approximants. Only close to Tq we find differences of the order of 10%. 

The expansion coefficients Cn,(T, m^) have been calculated in Q for m = 0.1 at a set of 16 temperature values in 
the interval T/Tq S [0.81,1.98]. From this set of expansion coefficients we have calculated the partial derivatives c/ 
at temperature T by averaging over finite difference approximants for left and right derivatives at T. With this we 
have constructed the expansion coefficients for the energy and entropy densities. 


= 3cn{T,mq) + c'^(T,mq) , 

Sn = {3.- n)cn{T,mq) + c'^{T,rnq) . (8) 

Results for the 2’^'^, 4*^ and 6*^ order expansion coefficients of energy and entropy densities obtained in the approxi¬ 
mation discussed above are shown in Fig.^ Also shown in this figure are the corresponding expansion coefficients for 
the pressure {pn = c„). For temperatures larger than T/Tq ~ 1.5 all expansion coefficients satisfy quite well the ideal 
gas relations, and ef'® = Ssf ®/2. As noted already in Ref. (2 results for the 2”"^ and 4*^ order expansion 

coefficients are close to those of an ideal Fermi gas which describes the high temperature limit for the energy density 
of 2-flavor QCD, 


10 V T / 27r2 V T / 


(9) 


Nonetheless, quantitative statements on the approach to ideal gas behavior may be at present premature. Deviations 
from the continuum result clearly arise in our current analysis which has been performed on lattices with a lattice 
spacing in units of the inverse temperature given by aT = 1/Nr = 1/4. Although cut-off effects are reduced due the 
usage of an improved action the approach to the continuum limit has to be analyzed further through calculations on 
lattices with larger temporal extent Nr ■ 

Although errors on the 6*^ order expansion coefficients are large it is apparent that subsequent orders in the 
expansion decrease in magnitude thus making the expansions for energy and entropy density well behaved for p,q/T<l. 
In fact, the 6*^ order contribution to e/T‘^ and s/T^ never exceeds 10% for Pq/T < 0.9. We will see below that this 
is the regime of interest for a discussion of the thermodynamics of matter formed in heavy ion collisions at the SPS 
and RHIC and even is appropriate for the energy range to be covered by the AGS and the future FAIR facility at 
GSI/Darmstadt. Results for the /Tg-dependent contribution to the energy and entropy densities, Ae/T^ = e/T^ — cq 
and As/T^ = s/T^ — sq, obtained in 4*^ order Taylor expansion are shown in Fig.[3for some values of Pq/T. 


III. LINES OF CONSTANT ENTROPY PER BARYON NUMBER 


When discussing the equation of state at finite baryon number density we have to eliminate the quark chemical 
potential in the thermodynamic quantities introduced in the previous section. We may do this by demanding that 
a thermodynamic quantity, for instance the quark number density, stays constant. We will use here a prescription 
more appropriate for discussing the thermodynamics of matter created in relativistic heavy ion collisions. After 
equilibration the dense medium created in such a collision will expand along lines of constant entropy per baryon. It 
then is of interest to calculate thermodynamic quantities along such isentropic lines. 

We thus should first analyze how the quark chemical potential changes along lines of constant S/Nb- It is readily 
seen that in an ideal quark-gluon gas these lines are defined by Pq/T = const, 


c 377r=^ I 

_ 2 45 V T / 

Nb /a j _If 

T ^ \ T ) 


( 10 ) 


In the zero temperature limit, however, the resonance gas reduces to a degenerate Fermi gas of nucleons and the 
chemical potential approaches a finite value to obtain finite baryon number and entropy, i.e. Pq/T ~ 1/T. 
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FIG. 1: The 2"''*, 4*^ and 6*^ order Taylor expansion coefficients for pressure {p„ = Cn), energy density (e„) and entropy density 
(sn) as functions of T/Tq. 



FIG. 2: The pq dependent contribution to energy density (filled symbols) and entropy density (open symbols) calculated in 
4*^ order Taylor expansion. Data points for the entropy density have been shifted slightly for better visibility. Horizontal lines 
show the corresponding ideal gas values for the energy density. 


In Fig. 0 we show the values for p,q needed to keep S/Nb fixed. These lines of constant S/Nb in the T-piq plane 
have been obtained by calculating the total entropy density for 2-flavor QCD from Eq. fusing results for the pressure 
and energy density calculated at /ig = 0 and the corresponding p,q dependent contributions shown in Fig. ^ 

The ratio of s/T^ and Uq/T^ obtained in this way is then solved numerically for piq/T. 

We find that isentropic expansion at high temperature indeed is well represented by lines of constant p.q/T down to 
temperatures close to the transition, T ~ 1.2ro. In the low temperature regime we observe a bending of the isentropic 
lines in accordance with the expected asymptotic low temperature behavior. The isentropic expansion lines for matter 
created at SPS correspond to S/Nb — 45 while the isentropes at RHIC correspond to S/Nb — 300. The energy range 
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FIG. 3: Lines of constant entropy per quark number versus fiq/T (left) and in physical units using To = 175 MeV to set the 
scales (right). In the left hand figure we show results obtained using a 4*^ (full symbols) and 6*^ (open symbols) order Taylor 
expansion of the pressure, respectively. Data points correspond to S/Nb = 300, 150, 90, 60, 45, 30 (from left to right). The 
vertical lines indicate the corresponding ideal gas results, /iq/T = 0.08, 0.16, 0.27, 0.41, 0.54 and 0.82 in decreasing order of 
values for S/Nb- For a detailed description of the right hand hgure see the discussion given in the text. 
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FIG. 4: The complete energy density (left) and pressure (right) evaluated on lines of constant S/Nb using Taylor expansions 
up to 6*^ order in the quark chemical potential. Shown are results for S/Nb = 30, 45, 300. 


of the AGS which also corresponds to an energy range relevant for future experiments at FAIR/Darmstadt is well 
described by S/Nb — 30. These lines are shown in Fig. Olright) together with points characterizing the chemical 
freeze-out of hadrons measured at AGS, SPS and RHIG energies. These points have been obtained by comparing 
experimental results for yields of various hadron species with hadron abundances in a resonance gas |l3L Il4| . The 
solid curve shows a phenomenological parametrization of these freeze-out data M- In general our findings for lines 
of constant S/Nb are in good agreement with phenomenological model calculations that are based on combinations 
of ideal gas and resonance gas equations of state at high and low temperature, respectively Em. 

Our current analysis yields stable results for lines of constant S/Nb also for temperatures T<Tq and values of the 
chemical potential /is ~ 400 MeV. This value of tb is larger than recent estimates for the location of the chiral 
critical point in 2-flavor E and (2-|-l)-flavor QCD. We thus may expect modifications to our current analysis at 
temperatures below Tq once we include higher orders in the Taylor expansion and/or perform this analysis at smaller 
values of the quark mass. 

We now can proceed and calculate energy density and pressure on lines of constant entropy per baryon number 
using our Taylor expansion results up to 0{p^). In Fig. we show both quantities for the parameters relevant for 
AGS (FAIR), SPS and RHIG energies. At high temperatures the relevant values of the quark chemical potential in 
an ideal quark gas are = 0.82, 0.54 and 0.08, respectively. However, as can be seen from Fig. |3at temperatures 
T<1.8To the required values for the chemical potentials are significantly smaller. In particular, for the AGS (FAIR) 
energies {S/Nb = 30) we find /r^/T = 0.77 at T ~ I.STq. The 6*^ order Taylor expansion thus is still well behaved 
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FIG. 5: Equation of state on lines of constant entropy per quark number versus T/To (left) and in physical units using 
To = 175 MeV to set the scale (right). The solid curve in the right hand figure is the parametrization of the high temperature 
part of the equation of state given in Eq. II3 


at these values of the quark chemical potential. The dependence of e and p on S/Nb cancels to a large extent in the 
ratio p/e, which is most relevant for the analysis of the hydrodynamic expansion of dense matter. This may be seen 
by considering the leading 0{p'//) correction, 

P _ 1 _ 1 eg - 3po / 

e 3 3 eg V 

In Fig.j^we showp/e as function of T/Tg (left) as well as function of the energy density (right). In the latter case we 
again used Tg = 175 MeV to set the scale. The small dependence on S/Nb visible in Fig.EIleft) gets further reduced 
when considering the ratio p/e at fixed energy density (Fig. ISDright)). The softest point of the equation of state is 
found to be (p/e)mm — 0.075 and within our current numerical accuracy it is independent of S/Nb- However, the 
analysis of the temperature regime T<Tq at present clearly suffers from poor statistics and still needs a more detailed 
analysis. We also find that the equation of state for Tq < T < 2Tg is well parametrized by 

1.2 

1 + 0.5e fm^GeV 

We note, however, that this phenomenological parametrization is not correct at asymptotically large temperatures as 
it is obvious that Eq.^lwould lead to corrections to the ideal gas result p/e = 1/3 that are proportional to while 
perturbative corrections vanish only logarithmically as function of T. 
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eg - 3po eg 
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IV. CONCLUSIONS 

We have determined Taylor expansion coefficients for the energy and entropy densities in 2-flavor QCD at non-zero 
quark chemical potential. At present these expansion coefficients have been determined from calculations performed 
at one value of the bare quark mass which gives a good description of QCD thermodynamics at high temperature. In 
this temperature regime lines of constant entropy per baryon are well described by lines of constant Pq/T. On these 
isentropic lines we have determined the equation of state and find that the ratio of pressure and energy density shows 
remarkably little dependence on the ratio S/Nb- 

The regime close to and below the transition temperature Tg is not yet well controlled in our current analysis. 
Smaller quark masses and the introduction of a non-vanishing strange quark mass, kept fixed in physical units, will 
be needed to explore this regime as well as the approach to the 3-flavor high temperature limit in more detail. 
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